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ABSTRACT 

We study three transient X-ray sources, that were bright in the central region of M31 
galaxy in the year 2000. Observations with Chandra and XMM-Newton allowed us for the 
first time in the history of X-ray astronomy, to build light curves of transient sources in 
M31 suitable for studying their variability on a time scale of months and, in some periods, 
weeks. The three sources demonstrate distinctly different types of X-ray variability and 
spectral evolution. XMMU J004234.1-t-411808 is most likely a black hole candidate based 
on the similarity of its X-ray light curve and spectra to typical transient low-mass X- 
ray binaries observed in our Galaxy. The outburst of CXO J004242. 0+411608 lasted 
longer than a year, which makes the source an unusual X-ray transient. The supersoft 
transient XMMU J004319. 4+411759 is probably a classical nova-like system containing 
a magnetized, rapidly-spinning white dwarf. We estimate a total rate of X-ray transient 
outbursts in the central bulge of M31 to be of the order ~ 10 per year. The rate of the 
hard X-ray transients (~ 5 year~^) in the central part of the Andromeda Galaxy appears 
to be comparable to that of the central part of our own Galaxy. 



1. INTRODUCTION 

There are about 200 known X-ray binaries in 
our Galaxy. These contain either a neutron star or 
a black hole ( |van Paradijs 1995| ). About one third 
of them are classified as transients (see I'ianaka 



Shibazaki 1996| for a review). Some of the X- 
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(M31), the closest spiral galaxy to our own (we 
assume ~760 kpc distance throughout this Letter 
- van der Bergh 2000), are of special interest, be- 
cause M31 is in many respects similar to the Milky 
Way. M31 was observed previously with Einstein 
and ROSAT, and hundreds of X-ray sources were 
detected (li'rinchieri & I'abbiano 19911; IPrimini 
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et al. 1993 ; Supper et al. 1997 ; Supper et al 
200 1| ). This Letter discusses transients discovered 
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matically, but irregularly, with peak luminosities 



above 10^^ erg/s and typical decay timesorder a 
month. This last group represent the majority of 
known Galactic black hole candidates, and pro- 
vides the best opportunity to study the processes 
near black hole in relative proximity to the Earth. 
Galactic transients have been therefore a focus of 
substantial observational and theoretical effort by 
X-ray astronomers. Much less, however, is known 
about X-ray transients in other galaxies. Previous 
observations with Einstein and ROSAT were too 
sparse in time for any systematic monitoring of 
even nearby galaxies. The new capabilities pro- 
vided by Chandra and XMM-Newton started to 
make the difference. 
X-ray observations of the Andromeda Galaxy 

^Los Alamos National Laboratory, Los Alamos, NM 87545 



in the M31 core in recent sensitive observations 
by Chandra and XMM-Newton. 



2. OBSERVATIONS AND DATA ANALYSIS 

We analyzed a series of Chandra data consisting 
of seven observations with Advanced CCD Imag- 
ing Spectrometer (ACIS) and seven High Reso- 
lution Camera (HRC) observations. We used all 
publicly available observations of M31 central re- 
gion covering the period from Oct 1999 to July 
2000 (Table |l|). A more detailed description of the 
observations can be found in [Garcia et al. 200U| 
(hereafter GOO) and |DiStefano et al. 200 1| . 



The data were processed using the CIAO 
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v2.1.2'^ threads. We generated X-ray images of 
the central region of M31 (Fig. |l|), and performed 
a wavelet deconvolution of the images to detect 
and localize point sources. Because of calibra- 
tion uncertainties below 0.3 keV and high back- 
ground above 7 keV, we excuded ACTS data out- 
side these energy limits. To estimate energy fluxes 
and spectra, we extracted counts within source 
circular regions of ~ 3 to 6" radii (depending 
on the distance of the source from the telescope 
axis). Background counts were extracted from 
the adjacent source-free regions. Spectra of the 
bright hard sources (XMMU J004234.1+ 411808 
and CXO J004242.0-F411608) were fitted with an 
absorbed power law model using XSPEC v. 11^. 
The ACTS count rates were converted into en- 
ergy fluxes using the analytical fits to the spec- 
tra for brighter sources or PIMMS^ with standard 
parameters otherwise. For the HRC the source 
fluxes were estimated by PIMMS, based on the 
observed count rates and the spectral parameters 
measured by ACIS and XMM-Newton. 

The central region of M31 was observed with 
XMM-Newton on June 25 and December 28, 2000 
(^hirey et al. 200T| , h ereafter SOI; Psborne et al 



2001|] hereafter OOl; rTrudolyubov et al. 2001 



hereafter TOl) clS cl part of Performance Verifica- 
tion (PV) program (Table |I|). In our analysis we 
used the data from two European Photon Imag- 
ing Camera (EPIC) instruments: the MOSl and 
M0S2 detectors ( [I'urner et al. 200 1|) . In all ob- 
servations the EPIC instruments operated in full 
window mode (30' FOV) with a medium optical 
blocking filter. 

We reduced EPIC data with the XMM-Newton 
Science Analysis System (SAS v5.01)^. To obtain 
the spectra of individual point sources, we used a 
source circle of ~ 20 — 30" and took backgrounds 
from adjacent source-free regions. We used data 
in the 0.3 — 7 keV energy range, for consistency 
with Chandra data and because of uncertainties 
in the calibration of the EPIC instruments. This 
is the default energy band for all fluxes and lumi- 
nosities presented below. 

The energy spectra of the two hard 
sources (XMMU J004234.1+411808 and 
CXO J004242.0+41608) were fitted to an ab- 

•^http:/ / asc.harvard.edu/ciao/ 

''http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/index, 
^http: / /heasarc. gsfc.nasa.gov/Tools/w3pimnis. html 
^http://xmm. vilspa.esa.es/user/ 



sorbed power law, and an absorbed black 
body radiation model was applied to fit the 
spectrum of the supersoft transient XMMU 
J004319.4+411759 (Table |). Since the flux from 
the latter was negligible above 1.5 keV, we consid- 
ered only its 0.3-1.5 keV spectrum. EPIC-MOSl 
and M0S2 data were fit simultaneously, but with 
independent normalizations. 

3. RESULTS AND DISCUSSION 

3.1. Hard X-ray Nova XMMU 
JOO4234. 1+411808 

A new bright source, XMMU J004234. 1+411808, 
was discovered on June 25, 2000 during the first 
XMM-Newton observation of M31 (SOI, OOl, 
TOl). Our analysis of the archival Chandra data 
revealed the presence of a bright X-ray source at 
the position consistent with XMMU J004234.1+ 
411808 as early as June 21, 2000, and again on 
July 2 and July 29 (Table p. Based on the 
Chandra aspect solution, which is currently lim- 
ited by systematics to ~ 0.6" accuracy, we mea- 
sure the position of XMMU J004234.1+411808 
to be a = 00^2^34.44^ 6 = 41°18'09.7" (2000 
equinox) . 

The X-ray light curve of XMMU 
J004234. 1+411808 shown in Fig. |a. Combining 
the results of Chandra and XMM observations, 
we place an upper limit on the effective duration 
of the outburst of ^40 days. We have no mea- 
surements of the outburst rise, but it was shorter 
than 20 days and followed by a slower decay. A 
fit to decaying part of the light curve with an 
exponential gives a characteristic e-folding time: 
Tfoid = 17 ± 2 days. The source flux declined 
from ~ 3.3 x 10^^^ erg s^^ cm~^ on June 21 to 
~ 3.4 X 10^^^ erg s~^ cm^^ on July 29, which cor- 
responds to the luminosity drop from ~ 2.3 x 10^^ 
erg s"^ to ~ 2.3 x 10^^ erg s~^ From the ACIS 
observations prior to the outburst, and assum- 
ing the same spectrum as measured by XMM, we 
compute a 2a upper limit to the quiescent lumi- 
nosity of ~ 1.5 X 10'^^ erg s~^, ~ 20 times lower 
than maximum measured luminosity. 

The energy spectrum of the source measured 
with ACIS-S on July 2 is similar to one detected 
with EPIC-MOS a week earlier (Fig. ^a) and both 
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may be fit to an absorbed power law with a ~2.1 
(see Table H for fit parameters) . 

The energy spectrum and X-ray luminosity 
of XMMU J004234.1+411808 correspond to the 
hard/low spectral state of Galactic black hole 
candidates ( [I'anaka fc Lewin 1995| ). Moreover, 
the X-ray evolution of the source is typical for 
well established class of X-ray novae with fast- 
rise-exponential-decay (FRED- type) lightcurves 
( |Chen, Shrader, fc Livio 1997| ) . For many of these 
objects, radial velocity measurements proved that 
binary systems harbor black holes. We propose 
that XMMU J004234.1+411808 is an X-ray nova 
and the first stellar mass black hole candidate in 
M31. 

3.2. Transient source CXO JOO4242.O+4II6O8 

The X-ray source CXO J004242.0-I-411608 was 
discovered on December 1999 during the first 
Chandra observation of M31 (GOO) and remained 
detectable throughout whole Y2K (OOl, TOl). 
The source was not detected with XMM on June 



29, 2 001 ( I Sniroy 200l| ), which poto a lower limit of 
)0- 



^400 days to the outburst duration. 

The source exhibited remarkable stability of the 
X-ray luminosity and spectral shape throughout 
the outburst (Fig. ^b, c and Table 0). The aver- 
age energy spectrum may be approximated with 
an absorbed power law with a ~ 1.5 — 1.6. We 
note that the ACIS spectra of the source tend to 
be generally flatter than EPIC-MOS spectra (Ta- 
ble H). Best-fit values of the absorbing column 
are also systematically higher for the ACIS. Both 
discrepancies may be possibly attributed to the 
presence of a 10 — 15% pile- up in the ACIS spec- 
tra. 

Long (^1 year) plateaus were detected in light 
curves of several black hole candidates, whereas 
shorter plateaus seem more typical for neutron 
star systems (Chen et al. 1997). The hard spec- 
trum and long plateau make this outburst remark- 
ably similar to the long-lasting hard transient 
GRS 1716- 249 ( (^unyaev et al. 1994| ; [Revnivtsev 
let al. 19981 



(SOI, 001). Our analysis of May 26 Chandra data 
revealed the presence of a bright X-ray source at 
the position a = 00^43^19. 76^ 5 = 41°17'57.6" 
(2000 eq., error circle ~ 0.6"), consistent with pre- 
viously reported XMM position (OOl). The source 
remained active during the June 21 Chandra ob- 
servation. An energy spectrum obtained with 
EPIC- MOSl detector on June 25 is shown in Fig- 
ure ^d. It may be satisfactorily fit to an absorbed 
black body radiation model with a characteristic 
temperature of ~ 60 eV and absorbing column of 
~ 1.4 X 10^1 cm-2. 

The spectrum of XMMU J004319.4+411759 is 
similar to the Galactic supersoft sources and may 
be interpreted result of thermonuclear burn- 
ing of the accreted matter on the surface of the 
white dwarf (see |Kahabka fc van den Heuvel 1997| 
for a review) . The source proved to be a pulsator 
with a ~ 865s period (for more detailed discussion 
see 001). The transient behavior of the pulsator 
hints that it may be a classical Nova in the su- 
persoft X-ray spectral phase, several tens of days 
after the peak of the outburst ( [Kahabka fc van 
den Heuvel IMTt OOl). 



3.3. Supersoft transient 
XMMU JOO4319. 4+411759 



3.4. On the rate of transient outbursts 

Given a number of transient sources detected 
with Chandra and XMM, one can estimate the to- 
tal rate of transient outbursts in the central bulge 
of M31. Assuming a luminosity of > 10^'' erg s~^ 
and an average source outburst duration to be 
more than a month, and taking into account the 
spacing of the Chandra and XMM observations, 
we estimate a rate of ~ 10 year~^. We note a 
large uncertainty of this estimate, due to the lim- 
ited statistics. About ~ 50% of a total amount of 
transients detected in M31 are hard X-ray sources 
corresponding to ~ 5 outbursts per year. It is in- 
teresting to compare their number with a number 
of X-ray transient detected in the corresponding 
region of the Galaxy.^ Based on regular monitor- 
ing of the Galactic Center with RXTE and Bep- 
poSAX, the rate of hard transient outbursts in 
this region can be estimated as ~ 4 year~^, which 
is comparable to our estimate for the central part 
of M31. This is firther evidence for the similarity 
between the stellar populations of the Galaxy and 
M31. 



The X-ray transient source XMMU J004319.4+ 
411759 was discovered on June 25, 2000 by XMM 

^Very high extinction of X-rays towards the center of our own Galaxy prevents effective detection of supersoft sources 
with current instrumentation. 
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Table 1 

Chandra AND XMM-Newton OBSERVATIONS OF THE CENTRAL REGION OF M31. 



Date 

(UT) 


TJD" 


Observatory 


Obs ID 


Instrument 


Exposure** 
(ksec) 


13/10/1999 


464 


Chandra 


303 


ACIS-I 


9.0 


11/12/1999 


523 


Chandra 


305 


ACIS-I 


4.1 


27/12/1999 


539 


Chandra 


306 


ACIS-I 


4.1 


29/01/2000 


572 


Chandra 


307 


ACIS-I 


4.1 


13/02/2000 


587 


Chandra 


270 


HRC-I 


1.4 


16/02/2000 


590 


Chandra 


308 


ACIS-I 


4.0 


08/03/2000 


611 


Chandra 


271 


HRC-I 


2.4 


26/05/2000 


690 


Chandra 


272 


HRC-I 


1.2 


01/06/2000 


696 


Chandra 


309 


ACIS-S 


5.0 


21/06/2000 


716 


Chandra 


273 


HRC-I 


1.2 


25/06/2000 


720 


XMM-Newton 


100 


EPIC 


34.8 


02/07/2000 


727 


Chandra 


310 


ACIS-S 


5.0 


29/07/2000 


754 


Chandra, 


311 


ACIS-I 


6.0 


28/12/2000 


906 


XMM-Newton 


193 


EPIC 


12.4 



" - TJD - truncated Julian Date, TJD=JD-2451000 



- total exposure time 



CHANDRA HRC-I 
JUNE 21, 2000 
M31 



XMMU J004234.1+411808 



XMMU J0431 9.4+41 1758 

* ■ 




.CXO J004242.0+41 1 608 



NUCLEUS ' 



5 arcmin 



Fig. 1. — The central region of M31 as it appears in a 1.2 ks Chandra HRC-I observation (June 21, 2000). The 
transient sources are marked with arrows. The position of nuclear source is shown for comparison. 
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Table 2 

Best-fit model parameters of the energy spectra of transient sources { Chandra- ACIS and 
XMM-Newton-MOS data). Parameter errors correspond to la level. 



Absorbed Power Law (0.3 — 7 keV energy range) 



Date 


Photon 


Nh 


Flux" 






Instrument 


(UT) 


Index 


(xl022 cm-2) 




(d.o.f.) 




CXO J004242. 0+411608 


13/10/1999 


1.52 + 0.10 


0.25 + 0.04 


8.63 + 0.27 


7.76 + 0.24 


58.2(44) 


ACIS-I 


11/12/1999 


1.35 + 0.19 


0.13 + 0.07 


8.12 + 0.53 


6.46 + 0.42 


34.9(42) 


ACIS-I 


27/12/1999 


1.60 + 0.21 


0.17 + 0.08 


6.83 + 0.43 


6.00 + 0.38 


48.6(44) 


ACIS-I 


29/01/2000 
16/02/2000 


1.84 + 0.22 
1.45 + 0.17 


U-<J' -0.18 

0.20 + 0.06 


6.38 + 0.36 
8.59 + 0.44 


7.03 + 0.40 
7.33 + 0.38 


42.8(53) 
52.3(62) 


ACIS-I 
ACIS-I 


01/06/2000 


1.25 + 0.14 


0.11 + 0.03 


8.40 + 0.36 


6.50 + 0.28 


45.6(46) 


ACIS-S 


25/06/2000 


1.61 + 0.05 


0.10 + 0.01 


8.05 + 0.16 


6.59 + 0.13 


212.2(230) 


EPIC-MOS 


02/07/2000 


1.30 + 0.12 


0.10 + 0.03 


8.08 + 0.35 


6.23 + 0.27 


38.8(23) 


ACIS-S 


28/12/2000 


1.54 + 0.11 


0.06 + 0.03 


5.77 + 0.25 


4.39 + 0.19 


52.3(44) 


EPIC-MOS 


XMMU J004234.1+411808 


25/06/2000 


2.13 + 0.07 


0.07 + 0.02 


2.29 + 0.07 


2.02 + 0.06 


79.9(96) 


EPIC-MOS 


02/07/2000 


1.96 + 0.54 


0.07(/ixed) 


1.44 + 0.24 


1.15 + 0.19 


4(5) 


ACIS-S 




Absorbed Black Body (0.3 — 1.5 keV energy rani 


ge) 




Date 


Tbb 


Nh 


Flux^ 


Li 




Instrument 


(UT) 


(eV) 


(XIO^^ 






(d.o.f.) 




XMMU J004234. 1+41 1808 


25/0G/2000 


Gl ± 2 


0.14 + 0.02 


2.85 ± 0.07 


1.29 ± 0.04 


53.8(48) 


EPIC^-MOS 



- absorbed model flux in the 0.3 — 7.0 keV energy range in units of 10~ erg s~ cm~ 

^ - estimated unabsorbed X-ray luminosity in the 0.3 — 7.0 keV energy range in units of 10^^ erg s 

- absorbed model flux in the 0.3 — 1.5 keV energy range in units of 10""^^ erg s"^ cm"^ 

- estimated unabsorbed X-ray luminosity in the 0.3 — 1.5 keV energy range in units of 10^^ erg s 
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Candidate BH X-ray nova 
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TJD=JD-2451000. 
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Fig. 2. — X-ray flux histories of three transient sources in the 0.3 — 7.0 keV energy range. Upper panel: hard X-ray Nova 
XMMU J04234.1-I-411808. Middle panel: transient supersoft source XMMU J004319.4+411758. Lower panel: Chandra 
transient CXO J004242. 0+411608. Filled circles, open circles and open squares denote ACIS, HRC-I and EPIC-MOS 
data respectively. The X-ray luminosities are in units of xlO^'' erg s~^. Upper limits correspond to a 2a level. 
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XMM- Newton 



EPIC-MOSl 



T 



T 1 1 1 — r 

Jun 25, 2000 



T 



1 1 1 1 — r 

Jun 25, 2000 

CXO J004242. 0+41 1608 

Candidate BH transient 

+ 



0.1 



XMMU J004234. 1+411808 
Candidate BH X-ray nova ^ 



0.01 - 



0.001 



0.1 




I I I I I I 1 1 

Dec 28, 2000 

CXO J004242. 0+41 1608 

Candidate BH transient 




H — I I I I I 



1 1 \ — h+ 

Jun 25, 2000 

XMMU J004319. 4+411758 

Supersoft transient 



0.01 




0.001 ^ c 



J I I I 




Energy, keV 



Energy, keV 



Fig. 3. — Representative count spectra of transient sources. EPIC-MOSl data, 0.3 — 7 kcV energy range. Spectral 
model fits are shown for comparison. Panel a: hard X-ray Nova XMMU J04234. 1+411808, Jun. 25, 2000 observation, 
fit by absorbed power law model (a = 2.13, Nh = 7 x 10^° cm'^). Panel b: Chandra transient CXO J004242.0+411608, 
Jun. 25, 2000 observation, fit by absorbed power law model (a — 1.62, Nh = 1 x lO^'^ cm~^). Panel c: The same 
source, Dec. 28, 2000 observation, fit by absorbed power law model (a = 1.54, Nh = 6 x 10^" cm~^). Panel d: transient 
supersoft source XMMU J004319. 4+411758, Jun. 25, 2000 observation, fit by absorbed black body radiation model 
(Tbb = 61 eV, Nh = 1.4 x lO^^ cm'^). 



